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ABSTRACT 

It is becoming clear that the atmospheres of the young, self-luminous extrasolar gi- 
ant planets imaged to date are dusty. Planets with dusty atmospheres may exhibit 
detectable amounts of linear polarization in the near-infrared, as has been observed 
from some field L dwarfs. The asymmetry required in the thermal radiation field to 
produce polarization may arise either from the rotation-induced oblateness or from 
surface inhomogeneities, such as partial cloudiness. While it is not possible at present 
to predict the extent to which atmospheric dynamics on a given planet may produce 
surface inhomogeneities substantial enough to produce net non-zero disk integrated 
polarization, the contribution of rotation-induced oblateness can be estimated. Using 
a self-consistent, spatially homogeneous atmospheric model and a multiple scattering 
polarization formalism for this class of exoplanets, we show that polarization on the 
order of 1% may arise due to the rotation-induced oblateness of the planets. The de- 
gree of polarization for cloudy planets should peak at the same wavelengths at which 
the planets are brightest in the near-infrared. The observed polarization may be even 
higher if surface inhomogeneities exist and play a significant role. Polarized radiation 
from self-luminous gas giant exoplanets, if detected, provides an additional tool to 
characterize these young planets and a new method to constrain their surface gravity 
and masses. 

Key words: polarization - scattering - planets and satellites: atmospheres - 
stars:atmosphere. 



1 INTRODUCTION 



Several young, self-luminous gas giant planets have been de- 
tected by di rect imaging ([Chauvin et al.|2004[ [Marois et al 



2008 2010 Lagrange et al. 2010 Lafreniere et al. 20101 



around nearby stars. These objects are now being charac- 
terized by photometry and even spectrosc opy ([Bowler et ah 
20T0{ [Patience et al.|[20Tol [Currie et a"L]|2011| [Barman et 
al. 20111 in an attempt to characterize their atmospheres 



and constrain the planetary masses. In the next few years 
many more such planets are almost certainly to be detected 
by ground-based adaptive optics coronagraphs, such as the 
P1640 coronagraph on Palomar, the Gemini Planet Imager, 
and SPHERE on the VLT ([Beichman et al. ||2010[ ). 



The characterization of the mass of a given directly im- 
aged planet can be problematical, since such planets typ- 
ically lie at large star-planet separations (tens of AU and 
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greater) and are thus not amenable to detection by radial- 
velocity methods. Instead masses must be estimated either 
by comparison of photometry and spectroscopy to plane- 
tary evolutionary and atmospheric models or by their grav- 



itational influence on other planets or disk (e.g. Kalas et 



al. 2005 Fabrycky fc Murray-Clay 20101. Model compar 



isons as a method for constraining mass can be ambiguous, 
however. Evolution models which predict luminosity as a 
function of age have yet to be fully tested in this mass range 
for young planets and at very young ages (< 100 Myr) the 
model luminosity can depend on the unknown initial condi- 
tions (j Marley et al. , 2007 | [Fortney et al., 2008) . The masses 



of the planets around HR 8799 estimated by cooling mod- 
els are apparently inconsistent with standard model spectra 



( [Bowler et al]|2010| [Barman et"aL][2011| [Currie et a"L][2011 ) 

and can lead to rapid orbital instabilities if circular, face- 
on orbits are assumed ( [Fabrycky fc Murray-Clay |[2010[ ). 
Finally the mass of the planetary mass companion to the 
brown dwarf 2M1207 b ( [Chauvin et al.|2004[ ) inferred from 
fitting of spectral models to observed near-infrared colors 
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is discrepant with the mass inferred from the companion's 



luminosity and the age of the primary ( Mohanty et al. 2007 1 



Discrepancies such as these may arise because young 
exoplanets exist in a gravity-effective temperature {g,Tcs) 
regime in which both the evolutionary and atmospheric 
models have yet to be validated. Fits of photometry and 
spectroscopy to predictions of atmosphere models depend 
upon the veracity of the models themselves, which — in the 
Tefi range of interest — in turn sensitively depend upon model 
cloud profiles, which are as yet uncertain. Extensive expe- 
rience with fitting models to brown dwarf spectra and pho- 



tometry ( Gushing et al. 2008 Stephens et al. 2009 1 re- 



veals that while effective temperature can be fairly tightly 
constrained, gravity determinations are usually less precise, 
uncertain in some cases by almost an order of magnitude 
in g. While there are low gravity spectral indicators rec- 



ognized from surveys of young objects (Cruz et al. 2009 



Kirkpatrick et al. 2006 1 these have yet to be calibrated 



by studies of binary objects which allow independent mea- 
sures of mass. Ideally for a single object with a given radius 
R, evolution model luminosity, which (for a known parallax) 
constrains R^T^g would be fully consistent with {g, T^a) con- 
straints from atmosphere model fitting. But, as noted above, 
this is often not in fact the the derived luminosity, 

mass, and radii of the companion to 2M1207 b as well as the 
HR 8799 planets are not fully internally self-consistent with 
standard evolution models. 

Given the likely future ubiquity of direct detections of 
young, hot Jupiters and the clear need for additional in- 
dependent methods to constrain planet properties, we have 
explored the utility of polarization as an additional method 
for characterizing self-luminous planets. 

Polarization of close-in giant exoplanets whose hot at- 
mosphere favours the presence of silicate condensates, is dis- 



cussed by Seager, Whitney & Sasselov ( 2000) ) and by Sen- 



gupta & Maiti (20061. While these authors considered the 



polarization of the combined light from an unresolved sys- 



tem of star and planet, Stam, Hovenier & Waters (20041 



presented the polarization of the reflected light of a resolved, 
directly-imaged Jupiter-like exoplanet. Since the polarized 
light of a close-in exoplanet is combined with the unpolar- 
ized continuum flux of the star which cannot be resolved, 
the amount of observable polarization in such case is ex- 
trmeley low - of the order of magnitude of planet-to-star 
flux ratio. Polarization measurements of directly-imaged ex- 
oplanets in reflected light is also challenging. The removal 
of scattered light from the primary star must be precise 
in both polarization channels so that the planet's intrin- 
sic polarization (which is a differential measurement) can 
be accurately determined. In any case no extrasolar planet 
has yet been imaged in scattered light, an accomplishment 



that will likely require a space-based coronagraph (e.g., Boc- 



caletti et al. 20111. Measuring polarization of thermally 
emitted radiation — as we propose here — is also difficult but 
does not require a planet to be close to the star (where the 
starlight suppression is most difficult) so that it is bright in 
reflected light. Furthermore extrasolar planets have already 
been imaged which raises the possibility of polarization ob- 
servations. 

It is clear from comparisons of model spectra to data 
that most of the exoplanets directly imaged to date have 
dusty atmospheres ( Marois et al.|[2008 Bowler et al.|[2010 



Lafreniere et al.~||2010| [Barman et al.|[2011| |Currie et al.| 



2011 



Skemer et al. 



2011 1. Glear atmospheres lacking dust 



grains can be polarized, but only at blue optical wavelengths 
where gaseous Rayleigh scattering is important ( jSengupta| 



& Marley 2009 \ . Since even the hottest young exoplanets 



will not emit significantly in the blue, grain scattering must 
be present for there to be measurable polarization in the 



near-infrared where warm giant planets are bright (Sen- 



gupta fc Marley^ ^20 10) . There are two temperature ranges 



within which we expect a gas giant exoplanet to possess sig- 
nificant atmospheric condensates. The first is L-dwarf like 
planets (roughly 1000 < Tctr < 2400 K) where iron and 
silicate grains condense in the observable atmosphere. The 
lower end of this range in the planetary mass regime is as 
yet uncertain. The second temperature range occurs in cool 
planets with atmospheric water clouds (Tcff < 400 K). There 
have yet been no confirmed detections of such planets. Here 
we will focus on the first category since such objects are 
brighter, more easily detectable, and the comparison to the 
field dwarfs is possible. 

Although survey sizes are fairly small, linear polariza- 
tion of field L dwarfs has been detected. IMenard et al. I 
( |2002 1 and Zapatero Osorio et al. ( 2005 I both report that a 
fraction of L dwarfs, particularly the later, dustier spectral 



types, are intrinsically polarized. Sengupta & Marley (20101 



find that the observed polarization can plausibly arise from 
emission of cloudy, oblate dwarfs, although to produce the 
required oblateness (20% or more) the dwarfs must have 
fairly low gravity for a field dwarf {g ~ 300 ms~^) and rapid 
rotation. The required rotation periods are brisk, as little 
as 2 hours or less, but are compatible with observed ro- 



tational velocities in at least some cases (see Sengupta & 



Marley (20101 for a discussion). Sengupta fc Marley (20101 



further find that the near-infrared polarization is greatest 
at Teft ~ 1600 K where their model condensate clouds are 
both optically thick and still prominent enough in the pho- 
tosphere to maximally affect the polarization. 

Surface inhomogeneities can also give rise to a net polar- 
ization ( Menard fc Delfosse p OO"?) and experience from the 
solar system confirms that irregularly spaced clouds are to 
be expected. Both Jupiter's and Saturn's thermal emission 
in the five-micron spectral window is strongly modulated 
by horizontally inhomogeneous cloud cover and it would 
not be surprising to find similar morphology in the atmo- 
spheres of exoplanets. In the presence of surface inhomo- 
geneity, the asymmetry that produces the net non-zero disk- 
integrated polarization would increase and hence a combi- 
nation of oblate photosphere and surface inhomogeneity can 
give rise to detectable levels of polarization. 

Exoplanets are even better candidates than L dwarfs to 
have an oblate shape and be polarized. With a lower mass 
and roughly the same radius as a brown dwarf (and thus a 
lower gravity), a rapidly rotating planet can be significantly 
oblate and consequently produce a polarization signal even 
without surface inhomogeneities. Here we explore the condi- 
tions under which the thermal emission from a warm, young 
exoplanet may be polarized and consider the scientific value 
of measuring exoplanet polarization. We first look at the 
issue of oblateness, then present a set of cloudy model at- 
mosphere calculations relevant to planets amenable to direct 
detection and discuss under which conditions their thermal 
emission may be polarized. Finally we discuss our findings 
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and explore how the characterization of an extrasolar planet 
may be enhanced by polarization observations. 



2 YOUNG GIANT EXOPLANETS 
2.1 Evolution 

When giant planets are young they are thermally expanded 
and boast larger radii and smaller gravity. |Fortney et al.| 
( |2008[ ) have computed evolution models for gas giant plan- 
ets with masses ranging from 1 to 10 Mj for ages exceeding 
10^ yr and we can use their results to predict the oblateness 
of thermally expanded young Jupiters with various rotation 
rates. Those authors modeled two types of evolution models. 
The first variety, termed 'hot starts', was most traditional 
and assumed the planets formed from hot, distended en- 
velopes of gas which rapidly collapsed. This calculation is 
comparable to that of most other workers in the field. They 
also presented calculations for planets fomed by the core ac- 



cretion planet formation process (see Lissauer & Stevenson 



(20071) which (depending on details of the assumed bound- 



ary condition for the accretion shock) produces planets that 
are initially much smaller and cooler than in the 'hot start' 
scenario. 

For the calculations here we choose to use the 'hot 
start' evolutionary calculation. We do this for several rea- 
sons. First, these models provide a reasonable upper limit 
to the radius at young ages and thus bound the problem. 
Second, at the large orbital separations that will, at least 
initially, be probed by ground based adaptive optics coro- 
nagraphic imaging, the core accretion mechanism may be 
inefficient at forming planets. Thus the gaseous collapse sce- 
nario may be more relevant choice. Finally the three planets 
observed around HR 8799 are all much brighter than pre- 



dicted by the Fortney et al. (20081 cold-start, but not the 



hot-start, cooling tracks. 

Figure 1 presents model evolution tracks for non- 



irradiated giant exoplanets from Fortney et al. (20081. On 



this figure planets age from the right to the left as effective 
temperature falls, the planets contract, and their surface 
gravity, g, increases. The dashed lines denote isochrones. 
This figure guides our selection of atmosphere models 
to evaluate for polarization studies. Groundbased corona- 
graphic searches for planets are expected to foucus on stars 



younger than about 200 Myr (e.g., McBride et al. 20111 



From the figure we see that at ages of 10 to 200 Myr we 
expect exoplanets with masses falling between 1 and 10 Mj 
to have g roughly in the range of 15 to 200 ms~^. 



2.2 Shape 

Both Jupiter and Saturn are oblate. The fractional differ- 
ence, f — 1 ~ Rp/Re, between their equatorial and polar 
radii, known as oblateness, are 0.065 and 0.11 respectively. 
The extent to which their equators bulge outwards depends 
on their surface gravity, g, and rotation rate, Q,, as well as 
their internal distribution of mass. The Darwin-Radau rela- 
tionship ( Barnes fc Fortney [2003 1 connects these quantities 
of interest: 
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Here K = I /{MR^), I is the moment of inertia of the spheri- 
cal configuration, and M and i?e are the mass and equatorial 
radii. 

The relationship for the oblateness / of a stable poly- 
tropic gas configuration under hydrostatic equilibrium is 



also derived by Chandrasekhar ( 1933 1 and can be written 

as 

M'^R^ 



3 g 



(2) 



where C is a constant whose value depends on the polytropic 
index. 

The above two relationships provide the same value 
of oblateness for any polytropic configuration. Equating 
Eq. (1) and Eq. (2), we obtain 



C = 



^) +5 



(3) 



Substituting the value of K for a polytrope of index n 
gives the value of the corresponding C. For example, K — 
0.4,0.261,0.205,0.155,0.0754 for n = 0,1,1.5,2,3 respec- 
tively. The corresponding values of C derived by |Chan-| 



drasekhar (1933 1 (p. 553, Table 1). are 1.875, 1,1399, 0.9669, 



0.8612, 0.7716 respectively 

The interiors of gas giant planets can be well approxi- 
mated as n = 1 polytropes. For the observed mass, equito- 
rial radii, and rotation rates of Jupiter and Saturn, expres- 
sion (2) predicts, with n = 1, an oblateness of 0.064 and 
0.11, in excellent agreement with the observed values. Fig- 
ure 2 presents the oblateness computed employing Eq. (2) 
as applied to 1 and 10 Mj planets at three different ages, 
10, 100, and 1,000 Myr using the [Fortney et al.| ([20081) hot- 



start cooling tracks. Also shown is the oblateness (0.44) at 
which a uniformly rotating n — 1.0 polytrope becomes un- 
stable ( James [1964 1. Clearly for rotation rates comparable 
to those seen among solar system planets we can expect a 
substantial degree (/ > 0.10) of rotational flattening. As 
gas giants age and contract the same rotation rate produces 
much less oblate planets. However for young, Jupiter mass 
planets rotation rates of 7 to 10 hours can easily produce 
/ ~ 0.2 even for planets as old as 100 Myr. More rapid ro- 
tation rates may produce even greater degrees of flattening. 
L dwarfs, with much higher surface gravity, must have even 
more rapid rotation rates to exhibit even modest flattening 
( ,Sengupta fc Marley|2010| . 



3 POLARIZATION OF YOUNG EXOPLANETS 

To explore the degree of polarization expected for various 
planet masses and ages we considered a selection of one- 
dimensional, plane-parallel, hydrostatic, non-gray, radiative- 
convective equilibrium atmosphere models with sixty verti- 
cal layers ([Ackerman & Marley |[2001[ [Marley et al. [ [2002 



[Freedman et al. [|2008| ) for specified effective temperatures, 
800 < Teff < 1200 K and surface gravities g = 30 and 
lOOmsec"^. We focus on this apparently limited parame- 
ter range since all gas giant exoplanets with masses below 
10 Mj will have cooled below 1200 K by an age of 30 Myr 



(see Figure 1 and also Fortney et al. ( 2008 1). The median age 
for nearby (< 75 pc) young stars that are likely targets for 
planet imaging surveys is 50 Myr ( McBride et al.|2011 1. For 
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our study we choose a lower limit of 800 K, well below the 
Tcft at which Sengupta & Marley (2010) predicted maximal 
polarization for field L dwarfs. At such temperatures dust 
clouds, if present globally across the disk, will lie at high 
optical depth and we expect produce a smaller polarization 
signal than the warmer objects. 

Indeed 800 K is well below the field dwarf L to T transi- 



cally the off diagonal terms of the scattering phase matrix 



are described by White ( 1979 1 and are very similar to the 



tion temperature of 1200 to 1400 K (Stephens et al. (20091 



and references therein) by which point most signs of clouds 
have departed. However there exists growing evidence that 
there is a gravity dependence to the effective temperature 
at which clouds are lost from the atmosphere and certainly 
the planets such as those orbiting HR 8799 are still dusty 
at effective temperatures near 1000 K ( [Bowler et al][2010l ). 
Observation of a polarization signal in a cooler exoplanet 
would provide powerful evidence for atmospheric dust. 

Some of the more massive young exoplanets (Af > 
8Mj) may have gravities in excess of our lOOmsec"^ up- 
per limit, but as we show below little oblateness-induced 
polarization is expected at high gravity in this TcB range 



(see also Sengupta & Marley (2010l). For example a surface 
gravity of <; = 100 m see" and Teg = 1000 K approximately 
describes an 8 Mj planet at an age of 100 Myr while values 
of 30msec~'^ and 800 K are expected for a 2Mj planet at 
an age of 60 Myr. We choose these values and a few others 
to illustrate the parameter space and the sensitivity of the 
results to variations in gravity and effective temperature. 
Each model includes atmospheric silicate and iron 



clouds computed with sedimentation efficiency (Ackerman 
fc Marley [2001 1 feed = 2. Preliminary studies by our group 
suggest that even dustier models with /sed ~ 1 might be 
necessary to reproduce the HR 8799 planets. However our 
previous work ( Sengupta fc Marley|2010" I has demonstrated 
that while /scd = 1 atmospheres do show greater polariza- 
tion than /sod ~ 2, the difference is slight when integrated 
over the disk. Other cloud modeling approaches are reviewed 
by |Helling et al.| ( |2008[ ). Some of these alternative cloud mod- 
eling formulations, such as those employed by Helling and 



collaborators (e.g., Helling & Woitke 2006 Helling et al. 
|2008[ ), predict a greater abundance of small particles high 
in the atmosphere than the Ackerman & Marley approach. 
Such a haze of small particles could potentially produce a 
larger polarization signal than we derive here. Polarization 
measurements may thus help provide insight into the verac- 
ity of various approaches. 

As in [Sengupta fc Marley| ( [2010[ ) we employ the gas 
and dust opacity, the temperature-pressure profile and the 
dust scattering asymmetry function averaged over each at- 
mospheric pressure level derived by the atmosphere code in 
a multiple scattering polarization code that solves the ra- 
diative transfer equations in vector form to calculate the 
two Stokes parameter I and Q in a locally plane-parallel 



medium (^Sengupta & Marley 20091. For each model layer 



we fit a Henyey-Greenstein phase function to the particle 
scattering phase curve predicted by a Mie scattering calcula- 
tion. A combined Henyey-Greenstein-Rayleigh phase matrix 



( Liu & Weng 2006 1 is then used to calculate the angular 



distribution of the photons before and after scattering. In 
the near-infrared the contribution of Rayleigh scattering by 
the gas to the overall scattering is negligible and the scat- 
tering is treated in the Henyey-Greenstein limit with the 
particle phase function computed from Mie theory. Specifi- 



pure Rayleigh case. For the diagonal elements the Henyey- 
Greenstein elements are used. In the limit of the scatter- 
ing asymmetry parameter approaching zero the matrix con- 
verges to the Rayleigh scattering limit. Finally, the angle 
dependent I and Q are integrated over the rotation-induced 
oblate disk of the object by using a spherical harmonic ex- 
pansion method and the degree of polarization is taken as 
the ratio of the disk integrated polarized flux {Fq) to the 
disk integrated total flux (Fj). The detailed formalisms as 
well as the numerical methods are provided in [Sen gupta fc| 
Marley [pOOgf. 



Figures 3 and 4 illustrate typical input properties of the 
models employed here. Figure 3 shows a model temperature- 
pressure profile along with iron and silicate condensate grain 
sizes as computed by our cloud model. Figure 4 shows the 
mean layer single scattering albedo, lJq, and scattering asym- 
metry parameter, cosO, as a function of wavelength near the 
peak opacity of the cloud. Within strong molecular bands 
the single scattering albedo approaches zero since gas ab- 
sorption dominates over the cloud opacity. Below the cloud 
base and far above the cloud both the albedo and asym- 
metry parameters are essentially zero in the near infrared 
as gaseous Rayleigh scattering makes little contribution to 
the opacity at those wavelengths. Note that for the com- 
puted particle sizes the cloud is strongly forward scattering 
at wavelengths where cloud opacity dominates molecular ab- 
sorption in agreement with a recent study by [de Kok et al.[ 
(20Tll). 



4 RESULTS AND DISCUSSIONS 

Figure 5 presents the computed thermal emission and po- 
larization spectra of an approximately 10 Myr old 2 Jupiter 
mass planet assuming rotation periods of 5 and 6 hrs. The 
striking dependence of polarization on the rotation rate 
arises from the sensitivity of oblateness to rotation period as 
seen in Figure 2. Generally speaking the degree of polariza- 
tion is highest at those wavelengths of low gaseous opacity 
where the cloud is visible while at other wavelengths, inside 
of atomic and molecular absorption bands, flux emerges from 
higher in the atmosphere and is less influenced by cloud scat- 
tering. While the degree of polarization peaks at the shortest 
wavelengths shown (from the infiuence of gaseous Rayleigh 
scattering), there is very little flux at optical wavelengths. 
However in the near-infrared, where windows in the molecu- 
lar opacity allow flux to emerge from within the clouds, the 
computed degree of polarization approaches 1%. In these 
spectral regions the planets will be bright, the contrast with 
the primary star favorable, and thus the polarization may 
be more easily detectable at this level. Beyond about 2.2 /xm 
thermal emission emerges from above the cloud tops and 
thus there is no signature of the scattering and the net po- 
larization is near zero. This pattern of polarization is diag- 
nostic of atmospheric clouds and is easily distinguished from 
other sources of polarization, for example a circumplanetary 
disk. 

Figures 6 and 7 show warmer model cases for the same 
gravity with similar behavior. These cases would apply to 
quite young planets at an age of less than ten million years. 
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but illustrate that the degree of polarization does not dra- 
matically increase at higher effective temperatures. Figure 
8 shows the variation with gravity. With a fixed rotation 
period of five hours, models with g of 56 and 100 ms^^ 
show very little polarization at any wavelength. These mod- 
els would correspond to approximately 4 to 6 Jupiter mass 
planets at ages greater than 10 million years, perhaps typical 
of the planet types that may be directly imaged. The sensi- 
tivity of polarization to gravity seen in this figure illustrates 
the promise of polarization, in the right circumstances, to 
provide a new constraint on exoplanet mass. 

Figure 9 generalizes these trends, showing the predicted 
polarization in I and J bands as a function of the rotational 
period Prot- For a fixed surface gravity and viewing angle, i, 
the degree of polarization does not vary substantially within 
the range of T^s between 800 and 1200 K. The polarization 
profiles in both bands increase with decreasing rotation pe- 
riod and the polarization is generally greater in J than in / 



some transition dwarfs (Artigau et al. 112009 Radigan et al 



2011 1. Cloud patchiness is also observed in images of thermal 



band. As is the case for brown dwarfs (Sengupta & Marley 



20101, for a given rotation period the polarization decreases 



with lower i. 

All of the cases shown in Figure 9 have an oblateness 
less than 0.44, the stability limit for an n = 1 polytrope. 
For g — 30ms~^ the stability limit is reached at a rotation 
period of about 4 hours, slightly less than the lower limit 
shown on the figure. Such short rotation periods may in 
fact be a natural consequence of giant planet formation in a 
circumstellar binary as the angular momentum of accreting 



gas naturally produces rapid rotation rates ( Ward & Canup 



2010) 



We conclude that a self-luminous gas giant planet-even 
with a homogeneous cloud distribution-will exhibit notable 
polarization (greater than a few tenths of percent) in the 
near infrared if the planet is (1) cloudy, (2) significantly 
oblate, and (3) viewed at a favorable geometry. An oblate 
shape is the easiest to obtain at low masses and modest ro- 
tation rates or higher masses and more rapid rotation rates. 
Higher effective temperatures, which would produce more 



dust higher in the atmosphere and more polarization (Sen- 
[gupta fc Marley|2010[ ), are generally excluded by the evolu- 
tion for ages greater than a few million years. More massive 
planets, which take longer to cool, have higher gravity and 
thus a smaller oblateness and less polarization (Figures 2 
& 9) for a given rotation rate. Given these considerations 
we believe the cases we have presented here are among the 
more favorable for homogenous cloud cover. While we have 
not considered every combination of parameters, the models 
presented here along with perturbations of those models we 
have also studied lead us to conclude that uniformly cloudy 
planets will not present polarization greater than a few per- 
cent and polarization is most likely to be found for young, 
low mass, rapidly rotating planets. 

However inhomogeneous cloud cover, which we have not 
modeled, may also affect the polarization spectra. Indeed 
an inhomogeneous distribution of atmospheric dust (e.g., 
Jupiter-like banding) would not be unexpected. Such band- 



ing may provide further asymmetry ( Menard fc Delfosse 
|2004[ ) and hence increase (or even decrease) the net non-zero 
polarization. A non-uniform cloud distribution may be the 
mechanism that underlies the L to T-type transition among 
brow n dwarfs ([Ackerman fc Marley ||2001| [Burgasser et al.| 
2002 Marley et al. 2010 1 and variability has been detected in 



emission from Jupiter and Saturn taken in the M-band (five- 



micron) spectral region (e.g. ,Wcstphal_1969 Westphal et al 



1974 Orton et al. 1996 Baines et al 



2005 1 , so patchiness 



may indeed be common. Polarization arising from patchy 
clouds still requires the presence of some clouds of course, 
thus any polarization detection provides information on the 
presence of condensates and-by extension-constrains the at- 
mospheric temperature. 



5 CONCLUSIONS 

The next decade is expected to witness the discovery of a 
great many self-luminous extrasolar giant planets (Beich-| 
man et al. 2010) . The masses, atmospheric composition 
and thermal structure of these planets will be characterized 
by photometry and spectroscopy. For some systems, other 
constraints, such as dynamical interactions with dust disks 
or potential instabilities arising from mutual gravitational 
interactions (e.g., |Fabrycky fc Murray-Clay [2010 1 may also 
contribute. Here we demonstrate that measurable linear po- 
larization in / or J bands reveals the presence of atmospheric 
condensates, thereby placing limits on atmospheric composi- 
tion and temperature. Polarization of thermal emission from 
a homogeneously cloudy planet is most favored for young, 
low mass, and rapidly rotating planets. A diagnostic charac- 
teristic of cloud-induced polarization is that the polarization 
peaks in the same spectral bandpasses as the fiux from the 
planet because photons are emerging from within the cloud 
itself as opposed to higher in the atmosphere (Figures 5 
through 8). 

Assuming that our atmospheric and condensate cloud 
models are reasonably accurate, we conclude that any mea- 
sured polarization greater than about 1% likely can be at- 
tributed to inhomogeneities in the global cloud deck. While 
we have not considered every possible model case, we find 
that our most favorable plausible cases do not produce no- 



tably greater polarization. Other cloud models (e.g.. Helling 
fc Woitke|2006 Helling et al.|2008 1 which incorporate more 



small particles high in the atmosphere may well produce a 
different result, thus polarization may help to distinguish 
such cases. For a fixed rottion period, the oblateness and 
thus polarization increases with decreasing surface gravity. 
In such situations polarization may provide a new constraint 
on gravity and mass. However for gravity in excess of about 
50 ms~^ and for T^s < 1200 K (corresponding to planet 
masses greater than about 4 Mj) we do not expect detectable 
amounts of polarization. Warmer and higher gravity field L 
dwarfs, can show measurable polarization since such cloud 
decks are higher in the atmosphere. For directly imaged 
exoplanets, however, we do not expect to encounter such 
high effective temperatures. For exoplanets with plausible 
Teft < 1200 K, Figure 9 shows that even if the rotation pe- 
riod is as rapid as 4.5 hrs. and the viewing angle is 90° at 
which the polarization is maximum, the percentage degree 
of polarization in thermal emission is not more that a few 
times of lO"'^. 

The aim of our study was to better understand the 
information conveyed by polarization about the properties 
of extrasolar giant planets directly imaged in their thermal 
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emission. We have found that in some cases polarization can 
provide additional constraints on planet mass, atmospheric 
structure and cloudiness. Combined with other constraints, 
polarization adds to our understanding, although there re- 
main ambiguities. A study of the polarization signature of 
partly cloudy planets would yield further insight into the 
value of polarization measurements for constraining extra- 
solar giant planet properties. 
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Figure 1. Evolution through time in T^g - g space of non- 
irradiated, metallicity [AI/H] = 0.0 giant planets of various 
masses. Solid lines are evolution tracks at various fixed masses 
ranging from 1 to 10 Mj. Dashed lines are isochrones for various 
fixed ages since an arbitrary 'hot start' initial model | |Fortney et| 
|al.|2008| . 



Figure 2. Rotationally induced oblateness as a function of ro- 
tational period for 1 and 10 Mj planets with three different ages 
(10 - solid, 100 - dashed, and 1000 Myr - dotted). Horizontal line 
is the stability limit for n = 1 polytropes assuming solid body 
rotation. More rapidly rotating planets would form a triaxial el- 
lipsoidal shape and eventually bifurcate. 
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Figure 3. Atmospheric temperature (T, lower scale) and cloud 
particle size, r^g, (upper scale) as a function of pressure, P, for 
one adopted model atmosphere with T^ff = 1000 K, g = 30ms~^, 
and /scd = 2. From Figure 2 this corresponds to approximately 
a 2.5 Mj planet at an age of a few million years. Particle sizes 
are shown for iron and silicate (forsterite) grains, the two most 
significant contributors to cloud opacity. Shown is the "effective 
radius" (r^ff) of the particles, a mean size precisely defined in 
Ackerman & Marley (2001). 
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Figure 4. Scattering properties as a function of wavelength. A, 
of a model layer near the 1 bar pressure level for the model at- 
mosphere shown in Figure 3. Shown are the layer single scatter- 
ing albedo, UJo, solid, and the layer asymmetry parameter, cos6, 
dashed. In strong molecular bands gaseous absorption dominates 
over scattering, thus lowering the mean layer albedo. 
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Figure 5. The emergent flux (A) and the disk-integrated degree 
of linear polarization P{%) (B) of non-irradiated exoplanets at 
different wavelengths at viewing angle i = 90° (equatorial view). 
In (B), the top solid line represents the polarization profile for a 
rotational period Prot = 5 hr while the bottom solid line repre- 
sents that for 6 hr. Note that while the polarization can be high 
at blue wavelengths there is very little flux there. 
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Figure 6. Same as figure 5 but with T^g = 1000 K . This is the 
result for the model characterized in Figures 3 and 4. 
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Figure 7. Same as figure 5 but with T^ff = 1200ii'. 



Figure 8. The emergent fiux (A) and the disk-integrated degree 
of linear polarization (B) of non-irradiated exoplanets for a fixed 
Teff = lOOOK but for different surface gravities. In (B), the solid 
lines from top to bottom represent the polarization profiles for 
surface gravity g = 30, 56 and 100 ms~^ respectively. The dif- 
ference in the emergent flux (A) for a fixed effective temperature 
but surface gravity varying over this range is not noticeable on 
this scale. 




Figure 9. Scattering polarization profiles of non-irradiated exo- 
planets with different rotational periods. The solid lines represent 
the percentage degree of linear polarization in J-band while the 
broken lines represent that in I-band. A variety of model cases are 
shown, all assume f^^^ = 2. Cases are shown for viewing angle 
i = 90° (equator view) and i = 60 and 45°. 



